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Abstract. We present ^^CO J ~ 1^0 and J ~ 2— >1 observations of the low metallicity 
(12 + log(0/H) = 7.74) Local Group dwarf irregular galaxy WLM made with the 15 m 
SEST and 14 m FCRAO telescopes. Despite the presence a number of HII regions, we 
find no CO emission. We obtain low upper limits on the integrated intensity {Ico ^ 0.18 
K km s^^ for CO (1^0)). The non-detection is consistent with the result of Taylor, 

o 



^ ' Kobulnicky & Skillman (1998), that dwarf galaxies below a metallicity of 7.9 are 
^ [ not detected in CO emission. WLM shows that this trend continues for low metallicity 
^ ' galaxies even as their metallicities approach 7.9. These results are consistent with the 



models of the metal poor ISM by Norman & Spaans (1997). By comparing our CO data 
' with observations of star formation in WLM, we find evidence for a high CO to H2 
conversion factor. 



Key words: ISM: molecules - Galaxies: individual: WLM - Galaxies: ISM - Radio lines: 
galaxies 



1. Introduction 



The WLM galaxy is a dwarf irregular that sits on the edge of the Local Group. Its distance 



is 0.9 Mpc (Minniti & Zulstra 199£), close enough for detailed studies of the stellar 



populations both from the ground (Minniti & Zulstra 1996) and from HST (Dolphin 
|2000D . As is typical for dwarf galaxies, the metallicity of WLM is low, 12 + log(0/II) = 



7.74 (Skillman, Terlevich & Melnick 1989 ), while there is a comparatively large amount 
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of atomic hydrogen {Mhi — 4.7 x 10^ Mq; Huchtmeier, Seiradakis & Materne 1981) 



Morris & Lo ( 1978] ) observed a single position at the center of the galaxy in the CO J 



1— i-0 line, obtaining a 3cr upper limit of 0.25 K. 

CO has been difficult to detect in dwarf galaxies, a fact which is now most often 
attributed to their low metal content (e.g. Taylor, Kobulnicky, & Skillman 199^ ). Taylor 



et al. (1998) found evidence for a dependency of the integrated CO J = 1^0 line intensity, 
IcOj upon metallicity, such that only dwarfs with metallicities greater than ^ 8 on the 
12 + log (0/H) scale were detected, even with extremely low upper limits obtained on 
the more metal poor galaxies. The metallicity of WLM places it between the detected 
galaxies and the more metal poor ones, making it a logical candidate for a new search 
for CO emission. The tremendous improvements in receiver technology since the 1970's 



allow for much lower upper limits than could be obtained by Morris Lo (1978) and for 
observations of multiple positions within the galaxy. Here we report on our search for 
CO emission in the WLM dwarf irregular galaxy. 

2. Observations and Data Reduction 
2.1. SEST 

WLM was observed using the 15 m Swedish ESO Submillimeter Telescope (SEST) at La 
Silla, Chile, between 1999 January 10 - 16. Observations were carried out simultaneously 
in the CO J = 1^0 line at 115 GHz and the J = 2^1 line at 230 GHz. The beamsize 
at 115 GHz is 45", and at 230 GHz is 23". Using the Low Resolution Spectrometers, 
velocity resolutions of 1.8 km s^^ and 0.9 km s~^ were obtained at 115 GHz and 230 
GHz, respectively. The receivers were tuned to the systemic velocity of the galaxy, -116 
kms^^. The system temperatures varied from 330 to 470 K at 115 GHz and from 250 to 
420 K at 230 GHz. The pointing was checked at the beginning of each observing session 
by observing the nearby SiO maser source R Aqr. 

Observations were made of seven positions within the galaxy, each in proximity to the 



groups of HH regions observed by Hodge & Miller (1995). In other Local Group dwarf 



irregulars, such as IC 10 and NGC 6822, CO emission has been found in the vicinity of 



HH regions (e.g. Wilson 1992; Wilson 1995) 



2.2. FCRAO 

Observations were also obtained with the Five College Radio Astronomy Observatory 
14 m telescope in New Salem, Massachusetts from 1-3 December 1999, 30 December 
1999 - 2 January 2000 and 10 - 14 February 2000. The SEQUOIA focal plane array 
receiver, consisting of 16 pixels operating in the 85 to 115 GHz band, was used to observe 
the CO J = 1^0 line. The beamsize of the telescope at 115 GHz is 44", and the spacing 
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between pixels in the array is 88". The Focal Plane Array Autocorrelation Spectrometers 
(FAAS) were used to obtain a bandwidth of 80 MHz with 313 kHz (0.8 kms^^) channels. 
Because of the southern declination of WLM the system temperatures were higher in 
the FCRAO data than in the SEST data, ranging from 600 to 1000 K. The pointing 
was checked every 2 hours each session by observing the SiO maser sources R Cas and 
T Ceph. 

The telescope was pointed so that the central 2 columns of 4 pixels in SEQUOIA 
fell within the optical galaxy. The outer two columns, although outside most of the 
optical light from WLM, were well within the bounds of the HI distribution (Huchtmeier, 



Seiradakis & Materne 1981). Pixel 11 in the array was aimed to point at one of the 



prominent groups of HII regions identified by Hodge & Miller (1995). The telescope was 
operated in a staring mode, where the position was not changed to fill in the gaps between 
the SEQUOIA pixels. 

The SEQUOIA observations are complementary with the SEST data because they 
are not biased to observing only in the vicinity of HII regions. While HII regions are 
an indication of recent massive star formation, and hence likely places to find molecular 
gas, CO emission can be found in areas with no corresponding Ha emission (e.g. IC 10, 



Becker |l990t NGC 4214, Walter et al. gOOOj). Figure 1 shows the location of the SEST 
and FCRAO observations compared to an optical image of WLM. 



2.3. Data Reduction 

The data were reduced using the CLASS package. We subtracted a linear baseline from 
each individual scan; scans for which a linear baseline was not sufficient were discarded. 
The final spectra for each position were obtained by averaging all the scans at each 
position. To decrease the rms noise in the final SEST spectra, we smoothed them to 
velocity resolutions of 7.2 and 16.2 kms^^ for the 1^0 data and to 2.7 and 7.3 kms^^ for 
the 2— >1 data. To convert to the main beam temperature scale, main beam efficiencies 
{tImb) of 0.70 for 115 GHz and 0.5 for 230 GHz (SEST Handbook) were used. The 
FCRAO spectra were smoothed to 7.2 kms~^, and rjMB = 0.45 was used. All the analysis 
was done using the spectra at 7.2 kms^^ resolution. 



3. Results 

Both the original resolution and the smoothed data were searched for CO emission. In 
neither the SEST data at either frequency, nor the FCRAO 115 GHz data was any 
emission detected. Table 1 gives the properties of the SEST data from the individual 
positions, including the offset in arcseconds from the center position, the rms noise in 
the spectra, and the 5cr upper limit on integrated intensity (Ico)j for both the (1— >0) 
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Fig. 1. Location of the SEST (hexagons) and FCRAO (circles) superposed over a digi- 
tized POSS image of WLM. 

and (2^1) data. Table 2 gives the rms noise and (Ico) 5cr upper limits for the FCRAO 
data. In all cases the upper limits are calculated assuming undetected emission at the 5a 
level in one 7.2 kms~^ channel. 
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Fig. 2. Example spectra from each telescope. The upper spectra show the original velocity 
resolution, while the lower are smoothed to 7.2 kms^^. The horizontal line shows the full 
width at 20% of maximum velocity range of the HI line from Huchtmeier et al. (|l98l| ). 
SEST-A and FCRAO-07 are the same location, as are SEST-E and FCRAO-11. 



4. Discussion 



To compare our WLM data with other gas rich dwarf galaxies, in Figure 3 we reproduce 
the plot of Ico versus galaxy metallicity from Taylor, Kobulnicky and Skillman ( 1998| ), 



with a data point for WLM added. The data point we added represents position A 
from the SEST observations. This was chosen because its spectrum is typical, being 
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Table 1: SEST Positions in WLM 



Position 


ARA 


ADec 


(TfT„h) Cl-O) 


Inn (1-0) 


a{Trr^b) (2-1) IcO (2-1) 






(arcsec' 


(arcsec) 


(K) 


(K kms"^) 


(K) 


(K kms"^) 


A 










0.005 


< 0.18 


0.003 


<0.11 


B 







-45 


0.004 


< 0.14 


0.003 


<0.11 


C 




-45 





0.004 


< 0.14 


0.003 


<0.11 


D 




-31 


-31 


0.004 


< 0.14 


0.003 


<0.11 


E 







-86 


0.009 


< 0.32 


0.005 


<0.18 


F 




-31 


-100 


0.004 


< 0.14 


0.003 


<0.11 


G 







-131 


0.004 


< 0.14 


0.002 


<0.07 






Note: The (0,0) position is 23 


'' 59"" 23.0^ - 


15° 43' 42.5". 








Table 2: FCRAO Positions in WLM 






Position G 


(T^t,) (1-0) 


Ico (1-0) 


Position o-(T^6) (1-0) 


Ico (1-0) 








(K) 


(K kms^^) 




(K) 


(K kms"^) 




1 




0.011 


< 0.40 


9 


0.016 


< 0.58 




2 




0.009 


< 0.32 


10 


0.013 


< 0.47 




3 




0.011 


< 0.40 


11 


0.009 


< 0.32 




4 




0.007 


< 0.25 


12 


0.016 


< 0.58 




5 




0.011 


< 0.40 


13 


0.011 


< 0.40 




6 




0.013 


< 0.47 


14 


0.013 


< 0.47 




7 




0.018 


< 0.65 


15 


0.009 


< 0.32 




8 




0.011 


< 0.40 


16 


0.007 


< 0.25 






Note: Position 11 corresponds to 23'' 59"" 23.0" 


, -15° 45' 


10.9". 



neither the lowest nor highest noise, and because this location coincides with large HII 
regions, and so is a likely place to have expected CO emission. Taylor et al. noticed in 
their data a trend for lower metallicity galaxies to be undetected while those of higher 
metallicity were detected. This was not a new result, but because they used a relatively 
large number of galaxies, they found that the division between detected and undetected 
galaxies occurred at 12 -f log(0/H) ~ 7.9. (about 1/10 solar abundance). WLM falls just 
below this metallicity at 7.74 and in a region in the figure with a dearth of data points. 
Our non-detection of WLM at the level of Ico < 0-18 K kms^^ in CO (1— >-0) reinforces 
this result. Our new data point shows that the trend for galaxies with metallcities below 
7.9 to be undetected in CO emission continues almost right up to a metallicity of 7.9. 



Norman & Spaans (1997) have studied the formation of molecular gas at high red- 
shifts, where the metallicity is low. They find that the cold, dense phase of the ISM does 
not form until a metallicity of approximately 0.03 to 0.1 Zq. In their models, before this 
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time star formation occurs at very low levels. At this metallicity C, O, and CO become 
abundant enough to contribute to cooling the gas, and giant molecular clouds can form. 
Because dust formation requires metals, it is also in this metallicity range that there 
is sufficient dust to shield molecular gas. Also, the amount of molecular gas increases 
because of the formation of H2 on the surface of dust grains. Although they made their 
models with protogalactic disks (presumably becoming spiral and elliptical galaxies with 
evolution) , the phase evolution of the ISM broadly fits what is known about dwarf galax- 
ies. Dwarf galaxies exhibit the lowest metallicities known from emission line regions, down 
to Z < 0.03 Zq (e.g. I Zw 18; Skillman & Kennicutt 1993| ). The lowest metallicities seen 
in these galaxies are close to the lower range of metallicity at which Norman & Spaans 
argue that GMCs can form. This would make sense if large star formation episodes can 
only occur in gas cool enough and dense enough that it is molecular. Our non-detection 
of WLM agrees with the Norman & Spaans models as well. The upper portion of the 
range 0.03 to 0.1 Z© corresponds to 12 + log(0/H) ~ 8.0, the metallicity at which dwarf 
galaxies are no longer detected in CO in Figure 3. The metallicity of WLM, 7.74, is below 
this level. 



Taylor et al. (1998) attributed the non-detection of dwarfs at metallicities below 12 
-I- log(0/H) ~ 8.0 to the CO to H2 conversion factor (Xco) having a dependence on 
metallicity which becomes sharply non-linear at this metallicity. However, there are not 
enough data points in Figure 3 to say if there is truly a sharp cutoff in the detectability, or 
if the small numbers make it appear so. Clearly more observations, especially of galaxies 
with metallicities in the range 7.8 < 12 -t- log(0/H) < 8.0, are needed before the behavior 
of CO emission at low metallicities can be understood. 



Verier & Hodge ( |1995|) obtained very low upper limits on CO emission in the nearby 
dwarf galaxy GR 8, and compared these limits to the H2 column density (Njjj) inferred 
from the presence of star formation in the form of HII regions. We will apply their analysis 



to WLM, using the Ha data of Hodge & Miller ( |1995| ): 

Nh. = 9.5 X 10-3 LMer9s'^)rsAyr) 

where Lua is the Ha luminosity, tsf is the timescale of star formation, esF is the star 
formation efficiency, and Rmc is the radius of a molecular cloud. For tsf we adopt 10^ 
yr, as did Verier & Hodge. This number ultimately derives from Lo, Sargent & Young 
(|199|), who argue that tsf must be between the cloud-cloud collision time scale (the 
timescale on which the gas will collapse to form stars) and the timescale for disruption of 
the ISM by OB stars. The former has an upper limit based upon the crossing time, which 
for dwarfs is ^ 2 x 10^ yr, while latter will take place on a time scale of a few x 10^ yr. 
Verier & Hodge assumed £5^ = 110 obtain lower limits on N//^ , but a more accurate 
value is 0.02 (e.g. Taylor et al. 1999] for the post-starburst dwarf galaxy NGC 1569, 
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Fig. 3. 

Taylor, 
added. 



Log(Ico) vs. oxygen abundance, 12 + log (0/H). The plot is reproduced from 



Kobulnicky & Skillman (1998), with the data point for WLM from this paper 



Wilson & Matthews |1995| for the giant HII regions NGC 604 and NGC 595 in M33). 
We will estimate the column density with both efficiencies since this is the least secure 
property used in the calculation. For Kmc we adopt 10 pc based upon spatially resolved 
observations of giant molecular clouds in nearby dwarf galaxies with interferometers (e.g. 



IC 10, Wilson & Reid |l99l| and Wilson |l995f NGC 6822, Wilson |199§. 

We can also estimate N^f^ from our CO data, if we assume a CO to H2 conversion 
factor. For this we apply the conversion factor - metallicity relation derived by Wilson 
( |1995D from Local Group galaxies, extrapolating down to the metallicity of WLM. This 
gives a conversion factor of 1.3 x 10^^ cm^^ (K kms^^)^^, or 5.8 times greater than 
the Galactic value. Note that if the conversion factor behaves as suggested by Taylor 



et al. (1998), then the true conversion factor could be much higher than this estimate. 
Therefore, our N/f^ estimates should be considered lower limits. 

In Table 3 we give the Ha luminosity, the H2 column density estimated from star 
formation (N/f.-, (SF)) for esF — 0.02 and 1.0, and the H2 column density lower limit from 
the CO observations (Nij2(C0)) for the four SEST positions which overlap with HII 
regions from Hodge & Miller ( 19951 ). For each position, the limits based on the CO data 
are already lower than the densities inferred from the star formation. This implies that 



the conversion factor is even higher than the extrapolation of the Wilson (1995) relation 
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would suggest. In order to get the average N//2(C0) up to the level of the average 
(SF) for the csf = 0.02 case would require a conversion factor more than 1500 times 
the Galactic value. Clearly this is extreme, and may result from assuming too low a value 
for 65^. In the case oi esp = 1 (which gives the minimum amount of molecular material 
for the observed star formation rate), a conversion factor of ^ 30 times Galactic value 
would be necessary. A more realistic esp would increase this number. Through the same 
analysis, Verter & Hodge arrive at a lower limit for the for the conversion factor in GR 8 
of ^ 30 times the Galactic value. GR 8 has a lower metallicity than WLM (12 + log 
(O/H) ~ 7.5 compared to 7.7), so in principle the conversion factor ought to be higher 
than in WLM, but the fact that we obtain the same upper limit here as Verter & Hodge 
merely represents that both galaxies have similar levels of star formation and we have 
reached similar sensitivities in the observations. 

Some dwarf galaxies have regions of CO emission with no corresponding HII regions, 



areas in a pre-star formation phase (Becker 1990, Walter et al. 2000). Such regions are 
not as common, however, as cases where the molecular gas and HII regions are in close 
association. Clearly in the case of a molecular cloud with no star formation, one could not 
expect to estimate the molecular gas content associated with newly formed stars. The 
crude calculation we describe above tells us nothing about such clouds. But HII regions 
in the Galaxy are often associated with molecular gas on large scales, e.g. the Orion 
region. This is also true of prominent HII regions in dwarf galaxies like IC 10 (Becker 
1990| ), NGC 1569 (Taylor et al. |l999D and NGC 4214 (Walter et al. |2000D . It is unlikely 



that all of the HII regions in WLM have totally disrupted their natal molecular gas, so 
it is reasonable to attempt to detect it, and to estimate the emission from these clouds 
based upon the young stars we currently see. 

Of course there are large uncertainties in the above calculation: in addition to the 
star formation efficiency, one must assume an IMF (Verter & Hodge use a Salpeter 
IMF), and neglect dust (observations of IC 10 have shown that dwarf irregulars can host 
considerable dust - Bolatto et al. |2000D . Substantial dust obscuration would mean we 
have underestimated the molecular gas implied by the observed star formation, which 
would further increase the conversion factor. Our data certainly do not prove a high 
CO-H2 conversion factor, but they are consistent with it. 



Klein (2000) has recently argued that the conversion factor, Xco, depends upon the 
cosmic ray flux. He used the 6 cm radio continuum surface brightness as a tracer of 
cosmic rays and found that as the cosmic ray content increased, Xco decreased. The 
implication here is that the heating by cosmic rays provides a significant portion of the 
energy powering the CO emission, so that more CO emission per unit H2 mass is observed 
when the cosmic ray flux is large compared to when it is small. WLM fits into this picture 
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Table 3: Star Formation Rates and Inferred H2 Column Densities 



Position L(Hq) Nff2(SF) Nh2(SF) Nn^iCO) 

esF = 0.02 esF = 1 

(10^^ erg s-i) (10^^ cm-2) (IQ^i cm-^) (IQ^i cm'^) 

A 2.0 95 1.9 > 0.24 

C 1.0 48 0.95 > 0.18 

E 0.8 38 0.76 > 0.42 

G 2.1 100 2.0 > 0.18 



well. It has a 20 cm continuum flux of only 5 mJy in the NVSS data (Condon et al. 199§| ), 
low compared to other dwarf galaxies with measured 20 cm fluxes (e.g. Hoppe et al. 1994, 
Klein, Weiland & Brinks 1991), suggesting a high Xco- The non-detection of CO is thus 
not unexpected. 



5. Summary and Conclusions 

We have searched for CO emission in the Local Group dwarf irregular galaxy WLM, 
and found none, obtaining low upper limits on the integrated intensity (Ico ^ 0-18 
K kms^^for CO (1^0). WLM has many HII regions, indicating recent star formation 
activity, so the lack of CO emission may indicate a high CO to II2 conversion factor, 
rather than a lack of molecular gas. The non-detection is consistent with the result of 



Taylor, Kobulnicky & Skillman ( |1998D that dwarf galaxies with metallicities less than 12 
+ log (0/H) ~ 8.0 are not detected in CO emission while those with greater metallicities 



are. This is explained in the models by Norman & Spaans ( 1997 ), who argue that CO 
is not abundant in the ISM until metallicities in the range ^ 0.03 to 0.1 Zq (12 + log 
(0/H) ~ 7.2 to 7.8). Comparing our upper limits to the molecular gas inferred from the 
presence of star formation in WLM suggests that the CO to II2 conversion factor may 
be as much as 30 times the Galactic value. 
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